Detection of the second eclipsing high mass X-ray binary in M 33 
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OO : ABSTRACT 

Chandra data of the X-ray source [PMH2004] 47 were obtained in the ACIS Survey of M 33 
(ChASeM33) in 2006. During one of the observations, the source varied from a high state to a low 
state and back, in two other observations it varied from a low state to respectively intermediate 
states. These transitions are interpreted as eclipse ingress and egresses of a compact object in a 
high mass X-ray binary system. The phase of mid eclipse is given by HJD 245 3997. 476±0. 006, 
the eclipse half angle is 30.6° ±1.2°. Adding XMM-Newton observations of [PMH2004] 47 in 2001 
we determine the binary period to be 1.732479±0. 000027 d. This period is also consistent with 
ROSAT HPJ observations of the source in 1994. No short term periodicity compatible with a 
rotation period of the compact object is detected. There are indications for a long term variability 
similar to that detected for Her X-l. During the high state the spectrum of the source is hard 
(power law spectrum with photon index ^0.85) with an unabsorbed luminosity of 2x 10 37 erg s _1 
(0.2-4.5 keV). We identify as an optical counterpart a V ~ 21.0 mag star with T e g > 19000 K, 
log(g) > 2.5. CFHT optical light curves for this star show an ellipsoidal variation with the same 
period as the X-ray light curve. The optical light curve together with the X-ray eclipse can be 
modeled by a compact object with a mass consistent with a neutron star or a black hole in a 
high mass X-ray binary. However, the hard power law X-ray spectrum favors a neutron star 
as the compact object in this second eclipsing X-ray binary in M 33. Assuming a neutron star 
with a canonical mass of 1.4 Mq and the best fit companion temperature of 33000 K, a system 
inclination i = 72° and a companion mass of 10.9 Mq are implied. 

Subject headings: Galaxies: individual: M 33 — X-rays: individuals: [PMH2004] 47 — X-rays: binaries 
— binaries: eclipsing 1 



1. Introduction 

Only a few eclipsing X-ray binaries (XRBs) 
have been detected in the Milky Way and other 
nearby galaxies. However, eclipsing XRBs are very 
interesting sources as they provide (when X-ray 
and optical data are included) the most reliable 
detailed mass measurements for compact stellar 
remnants possible (with the exception of the few 
known radio pulsars i n binaries, in cluding double 
neutron stars, see e.g. |Psaltidl2006ft . 

Of specific interest is the first eclipsing XRB 
in M 33, M 33 X-7 (hereafter X-7), which 
showed properties that were not known from 
XRBs in the Galaxy or the Magellanic Clouds. 
X— 7 was already detected as a variable by the 
Einstein observatory with a lu minosity at maxi- 
mum in excess of 10 38 erg s _1 ( Long et al.lll981 



Markert fe Rallisl Il983ft and it stayed active in 
all following observations. Its variability was ex- 
plained through the model of an eclipsing XRB 
with an orbital period of 3.45 d and an eclipse du- 
ration of ~0.4 d based on ROSAT and ASC A dat a 



Dubus et al 



(jLarson fc Schulmanl 11997 : 
With the help o f the well-constraine d Chandra 
position of X— 7, iPietsch et al.l (|2004aft identified 
a BOI to 071 star of 18.89 mag in V as optical 
counterpart: this star shows the ellipsoidal light 
curve of a high mass XRB (HMXB) with the X-7 
binary period. They argued that the compact ob- 
ject in the system is a black hole based on the 
mass of the compact object derived from orbital 
parameters and the optical companion mass, the 
lack of pulsations, and analysis of the extracted X- 
ray spectrum. Those authors concluded that X— 7 
would be the first eclipsing black hole HMXB. 

The b lack hole nature of X— 7 was firmly estab- 
lished bv lPietsch et al. ( 2006a ) using observations 
of the Chandra Advanced CCD Imaging Spec- 
trometer (ACIS) surv ey of M 33 (ChASeM33, see 
Plucinskv et al.ll2008 ) which accumulated in seven 
ACIS-I pointing directions a total exposure of 
200 ks each. During several of these pointings X— 7 
was in the field of view. The ChASeM33 measure- 
ments of X— 7 resolved for the first time the eclipse 
ingress and egress as well as constrained the light 
curve of X — 7 for binary phases around eclipse. 
In addition, IPietsch et al.l (|2006aft identified X— 7 
on archival HST WFPC2 images. Through de- 
tailed modeling of the optical light curve, the ob- 



served X-ray eclipse and improved parameters for 
the companion star, a lower mass limit for the 
compact object in the system of 9 Mq (see also 
Shporer et al .112007 ) has been obtained. The black 
hole mass was further constrained as 15.65 Mq 
with the help of a radial velocity curve based 
on Gemini North spectra and careful modeling 
( Orosz et al.1 12007). This - at the time - made 



X— 7 the stellar black hole with the highest mass 
that had been determined with high precision. 

In this paper we report on the detection and op- 
tical identification of the second eclipsing XRB in 
M 33 within the ChASeM33 project. The source 
was already det ected in the XMM-Ne wton sur- 
vey of M 33 of IPietsch et all (|2004bl hereafter 



PMH2004) and is source no. 47 in their cata- 
logue ([PMH2004] 47, hereafter PMH47). During 
the XMM-Newton observations it showed strong 
time variability and was classified as a tran sient 
XRB candidate by iMisanovic et all (|2006ft . It 
was first reported in the M 33 ROSAT cata- 
logue by lHaberl fc Pietschl (|200lft as source 19. 
PMH47 is listed as n o. 8 in the "first l ook" 
ChASeM33 catalogue (jPlucinskv et al.ll2008ft . In 
several ChASeM33 observations the source was 
detected far off axis with strong time variabil- 
ity which - in a preliminary analysis - could be 
explai ned as the signatur e of a 1.7 d eclipsing 



XRB (|Pietsch et alJl2006bft . A massive star from 



the Local Group su rvey catalogue of M 33 stars 
( Massev et al. 2006ft could be identified as the op- 
tical counterpart of the X-ray source. In a re- 
analysis of the deep 3.6m Canada-France-Hawaii 
Telesc ope (CFHT) photome tric survey images of 
M 33 ( Hartman et al. 2006ft . the proposed opti- 
cal companion showe d the ellipsoidal ligh t curve 
of a HMXB in M 33 (jShporer et al.ll2006ft . A co- 
herent analysis of the Chandra and XMM-Newton 
X-ray and optical data allows us to further con- 
strain the system parameters. Our follow-up anal- 
ysis of PMH47 in the ROSAT data also indicates 
time variability consistent with the Chandra and 
XMM-Newton ephemeris. 

2. Observations and data analysis 

The PMH47 area was covered by the ACIS 
CCDs during seven observations of the ChASeM33 
project and during the archival observation with 
observation identification (ObsID) 1730. Ta- 
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ble Q] summarizes these observations giving the 
ChASeM33 held in column 1, ObsID (2), obser- 
vation start date (3), sum of good time intervals 
("on-time", 4), the ACIS chip covering PMH47 
(5), the offset of the source from the pointing 
direction (6), and the PMH47 binary phase and 
cycle number during the observation (7, 8) using 
ephemeris that will be discussed in Sect.El In the 
last column (9), we comment on source brightness 
and observed binary features (see Sect. 2] and 0. 
The average brightness of PMH47 varied strongly 
between the ChASeM33 observations. However, 
the determination of accurate count rates and 
light curves turned out to be quite difficult be- 
cause in all observations the source is positioned 
far off axis and therefore has a large extent due to 
broadening of the telescope point spread function 
(PSF) which may stretch over bad CCD columns 
(e.g. in ObsID 6377) or may partly fall outside 
the detector field of view (FoV) (6382, 6385, 6387, 
7226, 7344). An extreme case is presented in Ob- 
sID 6383 where the position of the source falls 
outside the detector and only a small part of the 
PSF is recorded. In several of the observations, 
counts from PMH47 are rejected in standard level 
2 event files when the source is moving across re- 
jected columns due to satellite dithering. This 
effect can reduce the counts in 1000 s integration 
intervals by varying amounts and create spurious 
periods in timing analysis. Therefore, we created 
new level 2 event files for broadband time variabil- 
ity analysis that did not reject these columns. For 
better comparison we applied an approximate flux 
correction to account for the size of the extraction 
region relative to the PSF size and to convert to 
an equivalent on-axis so urce. For each obs erva- 
tion, a CriaRlQ raytrace ( Carter et al. 2003) was 



performed for a 1.5 keV on-axis point source for 
an input ray density of 1 mm -2 , and the total 
number of resulting rays was evaluated. For each 
observation, another ChaRT run was performed for 
a 1.5 keV point source (ray density also 1 mm~ 2 ) 
at the appropriate off-axis position; this applies 
the mirror reflectivity and vignetting. The rays 
were projected to the detector using the CIAO tool 
psf .project jray, which also applied an approxi- 
mate detector quantum efficiency (QE). In many 
cases, the source is near the edge of the detector, 



or in one case, off the detector (although some of 
the large PSF at that position is still on the de- 
tector). To account for these exposure variations 
resulting from chip edges and spacecraft dither, we 
weighted "counts" images generated from the ray- 
trace pseudo-event lists by a fractional exposure 
map. The fractional exposure map accounts for 
the fraction of the time a sky pixel is on the detec- 
tor, and is obtained by constructing an "exposure 
map" in which the Chandra High Resolution Mir- 
ror Assembly (HRMA) effective area and detector 
QE are both set to unity. We extracted weighted 
counts from a region corresponding to the ex- 
traction region in the real observation to get the 
predicted number of counts within the extraction 
region used for the real data. This provides the 
predicted number of counts for the given source in- 
tensity. The normalization to an on-axis source is 
obtained by dividing the on-axis predicted counts 
by the extracted off-axis predicted counts. This 
factor was used to scale the observed fluxes to the 
values they would have had on-axis. For spec- 
tral fitting, PMH47 photons were extracted from 
CIAO v4.0 level 2 event files using the same area 
as for the time variability analysis, but with the 
bad columns excluded. Spectra were grouped in 
channels with at least 20 photons in the on-source 
spectra. 

In the XMM-Newton EPIC M 33 raster project, 
the PMH47 area was covered in seven observations 
at different off-axis angles (typical ontime of 14 
ks) spanning about three years. Table [2] summa- 
rizes these observations giving observation iden- 
tification (ObsID) in column 1, observation start 
date (2), on-time (3), the offset of the source from 
the pointing direction (4), 0.2-4. 5 keV luminosity 



(Col. 5, derived from Table 5 of iMisanovic et al 
2006), and the PMH47 binary phase and cy- 



1 http://cxc.harvard.edu/chart/ 



cle number during the observation (6, 7) using 
ephemeris that will be discussed in Sect.El In the 
last column (8), we comment on source bright- 
ness and observed binary features. The X-ray 
source was detected only during two observations 
in July 2001 which were separated by about 2.5 
days. However, during an observation between 
these two, the source was not detected. For the 
XMM-Newton EPIC X-ray spectral analysis we 
used PN (single + double pixel events, PATTERN 
0-4) and MOS (PATTERN 0-12) events disre- 
garding bad CCD pixels and columns (FLAG 0). 
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The EPIC spectra were simultaneously fit for the 
three instruments allowing for a constant normal- 
ization factor between the spectra. Spectra were 
grouped in channels with at least 30 photons in 
the on-source spectra. 



ing systematics ( Misanovic et al.l 20061) . In the 



ChA SeM33 "first look" paper (jPlucinskv et al 



In the ROSAT catalogue of M 33 (jHaberl fc Pietsch 



20011 ). the detection of PMH47 was reported in the 
HRI at a level of significance greater than la while 
in the PSPC it is only marginally detected at a 
level of significance of only 3<r. We reanalyzed the 
ROSAT HRI data of ObsIDs 600487h, 600488h 
and 600489h performed from July 27 to August 8, 
1994, where PMH47 is in the FoV. After screening 
for high background, we combined data with con- 
tinuous observation intervals that led to variable 
integration times of typically 1700 s (minimum 
786 s, maximum 2514 s), depending on the dura- 
tion of the scheduled observation and background. 

The PMH47 field was covered by the CFHT 
variability survey of M 33 conducted in the g', 
r', and i' bands during 27 nights from 2003 Au- 
gust to 2005 January. For details of this sur- 
vey and the d a ta red uction procedure used see 



Hartman et al 



(2006). To summarize the pro- 
cedure, the standard CCD calibrations were ap- 
plied to the images as part of the CFHT Queue 
Service Observing mode. We performed image 
subtraction photometry (jAlard k. Lupto n 1998; 
Alardl 2000l ) on these images to obtain differen- 
tial flux light curves for the source and PSF fit- 
ting photometry on a stacked reference i mage us- 
ing the DA OPHOT/ALLSTAR package (|Stetson 
1987lll992h to obtain the reference flux values. We 
then applied a scaling factor to the formal photo- 
metric uncertainties determined by requiring that 
X 2 per degree of freedom scatters about unity for 
~ 1000 non-variable stars spanning the full dy- 
namic range of magnitudes. 

For the data analysis we in addition used 
tools in the ESO-MIDAS v05SEPpll.0, EXSAS 
v03OCT_EXP, CIAO v3.2 and v4.0, and HEA- 
soft v6.3 software packages as well as the imaging 
application DS 9 v4.13 and the s oftware package 
ACIS Extract (iBroos et alj|2002f >. 



3. Improved position 

The best XMM-Newton position for PMH47 
was derived as RA J20 oo = 01 h 32 m 36 i !85, Dec J20 oo 
= +30°32'28'/5 with a la error of 0.57" includ- 



120081 ) the position of PMH47 is given as RAj2ooo — 
0l' l 32 m 36!95, Dec J2 ooo = +30 o 32'30'. / l with a la 
statistical error of 0.48". Assuming for PMH47 
the same systematic error as for sources within 3' 
from the aim poin10 a 99% error of 0.8" has to 
be added. We obtained the most reliable source 
position for ObsID 6385, as the vast majority of 
the PSF was contained on the ACIS-I array in 
this observation and the source was bright dur- 
ing most of the integration. Our ACIS Extract 
analysis of this observation determined the mean 
position of the source counts within the 90% en- 
circled energy PSF for the sources in the image, 
as well as the source position that provided the 
best match to a cross-correlation between the dis- 
tribution of source counts and the PSF shape at 
this location on the ACIS-I array. These two 
position determinations were in very good agree- 
ment for X— 7 and PMH47, both of which were 
contained in this observation. We then deter- 
mined the offset between the X— 7 position and 



its known optical counterpart in the IMassev et al 
(2006) catalogue. This offset we assume to be 



the systematic offset between the Chandra X-ray 
image and the IMassev et al. I (120061) catalogue im- 
ages. The offset was +0.5" in RA and +0.3" 
in Dec, with an alignment error of ±0.3". Af- 
ter making this correction, the final position and 
error determined for PMH47 in ObsID 6385 was 
RA,,2ooo = 01 h 32™36!94, Dec J2000 = +30°32'28"4 
with a la statistical error of ±0.37" correspond- 
ing to a la total error of 0.48". In Fig. 2]we show 
the improved po sition together with t he XMM- 
Newton position (jMisanovic et al. 2006) and that 
from the ChASeM33 "first look" paper overlaid 
on an optical imag e of the Local Group Survey 
(jMassev et al.ll2006l ). 



4. X-ray time variability 

For the ChASeM33 data, we sampled back- 
ground and solar system barycentre corrected light 
curves of PMH47 with a time resolution from 
1000 s to 5000 s depending on the brightness of 
the source in the observation. To increase the sig- 
nal to noise in the far off-axis observations, we re- 
stricted the analysis to the 0.5-5 keV band which 



2 see http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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covers most of the source flux (see Fig. [T]). In ad- 
dition we indicated in the images to the right of 
the light curves the extraction regions for source 
and background. Rectangular regions were used 
in some cases to get similar background structure 
in source and background regions. The dashed el- 
lipses show the HRMA PSF 90% enclosed counts 
fraction. For observations far off axis, we used ex- 
traction areas smaller than 90% PSF to improve 
the signal to noise ratio (e.g. for ObsIDs 1730). 
The light curve information has been confirmed 
by investigating time-selected images. In the fol- 
lowing we give a short description of the intensity 
behavior of the source within the individual obser- 
vations. 

During ObsIDs 1730 and 7344 no flux is de- 
tected from the source (3tr ACIS-I on-axis upper 
limit of 8xl0 -3 ct s _1 in the 0.5-5 keV band in 
5000 s integration intervals). 

During ObsID 6382, the source is faint and 
we needed to integrate the data over 5000 s. 
During the first 10 ks of the observation (untill 
HJD 2453698. 164±0. 029) no significant flux is de- 
tected. 

During the short ObsID 7226, the source is 
clearly detected but faint. The light curve again 
had to be integrated over 5000 s. No significant 
variability is detected. 

During ObsID 6383 flux at the wing of the 
PSF of the source is clearly seen in the im- 
age. However, from HJD 2453902.052±0.029 to 
HJD 2453902.342±0.029 no photons at all are de- 
tected at that position. 

During ObsID 6387, PMH47 is at least partly 
detected in a high state which allowed us to sam- 
ple the light curve with 1000 s resolution. We ob- 
served a transition of the source from zero to high 
intensity at HJD 2453912. 740±0.012 about 15 ks 
after the start of the observation. 

During ObsID 6385 PMH47 is again for most 
of the time detected at high intensity suitable 
for 1000 s binning. We observed transitions 
of PMH47 from high intensity to zero flux and 
back. The source transits to the low state at 
HJD 2453997.329±0.006 and returns to the high 
state at HJD 2453997.624±0.006. Ingress into and 
egress out of the low state can not be resolved by 
the 1000 s time bins imposed by source statistics. 
The duration of the low state is (25500±1000) s. 



After the low state the source seems to be more 
variable than before. Figure [2] demonstrates that 
during the low state no counts are detected from 
PMH47. 

During ObsID 6377, the source is not detected 
in the first -25 ks (until HJD 2454004.555±0.020) 
and then stayed in a faint state with indications for 
variability. Due to the faintness of the source the 
light curve was integrated with 3000 s resolution. 

In the bright phase of Chandra ObsIDs 6385 
and 6387 and in XMM-Newton ObsID 0102641101 
we searched for pulsations in the frequency range 
10~ 4 Hz to 7 Hz and found no significant peri- 
odic signal. For the Chandra observations, the 
2cr upper limit for a pulsation amplitude is 25% 
for frequencies below 0.15 Hz while for the XMM- 
Newton observation, the EPIC PN upper limit is 
43% for frequencies below 7 Hz. As one expects 
a modulation of the photon arrival time with the 
orbital phase which would smear out the pulsa- 
tion signal for short periods (amplitude of less 
tha n 50 s for HMXBs with short orbital period, 
c.f. iBildsten et al.lll997l ). we not only investigated 
power spectra for the entire observation, but also 
added up power spectra for intervals as short as 
3319 s (1024 time bins with the Chandra instru- 
ment resolution of 3.241 s). 

5. Energy spectra 

For the brighter states of PMH47 that is, dur- 
ing the high states of Chandra ObsIDs 6385 and 
6387 as well as XMM-Newton ObsID 0102641101 
- more than 900 counts were collected from the 
source, allowing us to perform a detailed spec- 
tral analysis. The resulting integration times, raw 
count rates and degrees of freedom are listed in 
Table [3] together with the results of the power law 
fits. For XMM-Newton, luminosity refers to the 
EPIC PN spectrum, MOS2 yields similar values 
while those for MOS1 are ~7% higher. 

Due to the limited number of photons we only 
fitted absorbed one component spectral mod- 
els. We used two absorption components, ac- 
counting for the Galactic foreground absorption 
(with a fixed hydrogen column densi ty of 6xl0 20 
cm -2 and elemental abundances fromlWilms et al 



2000) and the M 33 absorption (with column den- 
sity as a free parameter in the fit and Galactic 
metal abundances). As expected for an X-ray 
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binary, thin thermal plasma spectra with abun- 
dances fixed to solar do not give acceptable fits. 
Also thermal bremsstrahlung fits can be rejected 
as the best fit temperature in all fits is only con- 
strained by the upper boundary (200 keV) mim- 
icking a hard power law. A similar argument holds 
for disk black body models. There, the best inner 
disk temperature values T- ln are around 5 keV and 
the lu minosity a few times 10 37 erg s _1 . Accord 



ing to lMakishima et aT ( 200df ). such temperatures 
are not reached even in ultra-luminous compact 
X-ray sources. Stellar mass black hole systems of 
similar luminosity to PMH47 show inner disk tem- 
peratures well below 1 keV. Power law fits for all 
observations have a similar or slightly worse mini- 
mum reduced \ 2 compared to the disk black body 
models. They are consistent with a common pho- 
ton index of 0.85. The Ng values for the power law 
model fits indicate some absorption within M 33 or 
intrinsic to the source in addition to the Galactic 
value (5.97 and 6.32 xlO 20 H cm" 2 in the direc- 
tion of PMH47 accordi ng to Dickey fc Lockmanl 
19901 : IStark et al.l ll992L respectively). Figure [3] 
shows the spectrum of the bright state of Chandra 
ObsID 6387 as an example. Un-absorbed source 
fluxes in the 0.2-4.5 keV band are in the range 
(2.4-2.7) xl0~ 13 erg cm~ 2 s _1 based on the best 
fitting power law model. These fluxes correspond 
to source luminosities of (1.8-2.0) x 10 37 erg s _1 , 
respectively. Throughout t he paper, we assume a 
distance to M 33 of 795 kpc ( van den Bergh|[l991 ). 

The JVh values of the power law model fits indi- 
cates that PMH47 lies in the plane of the M 33 disk 
or even on the near side, as the absorbing column 
within M 33 can be determined to 1.4-1. 9x 10 21 
cm~ 2 fr om a 47" x 93 " half power beam width 
Hi map (Newton 1980). Assuming total Ng val- 
ues from 0.6-1.5 xlO 21 cm -2 we can compute the 
expected optical extinction Ay to 0.34-0.84 mag 
and EjB - V) to 0.12-0 28 us ing standard rela- 
tions (jPredehl fc Schmittlll995h . 



6. Orbital period determination 

The intensity changes of PMH47 from high to 
low and vice versa in the Chandra ObsIDs 6385, 
6387 and 6377 are not resolved by the time reso- 
lution of our plots, of 1000 s, 1000 s, and 3000 s, 
respectively. Shorter time bins do not help as the 
source is not bright enough. In the following dis- 



cussion we try to interpret the intensity transitions 
reported in Sect. [4] by eclipses of a compact X-ray 
source by a companion star in a binary system 
and to determine eclipse parameters and orbital 
period of the system. For this analysis we use Ob- 
sIDs 6385, 6387 and 6377 and compare the results 
to the count rate behavior of PMH47 in the other 
observations. 

From observation 6385 we derive the eclipse du- 
ration as (25500±1000) s and a well defined epoch 
of mid eclipse (HJD 245 3997.476±0.006). How- 
ever, the orbital period is too long to detect a sec- 
ond eclipse during this 91 ks observation. 

A minimum length for the binary period can be 
estimated as 1.04 days by adding up the eclipse 
duration in ObsID 6385 and the longest continu- 
ous time out of eclipse (65000 s in ObsID 6387). 
To further constrain the orbital period of the sys- 
tem we use times of eclipse egresses. In HMXB 
systems the time of eclipse egress in most cases is 
better determined than the time of ingress which 
can be masked by additional absorption due to 
the viewing geometry through the innermost re- 
gions of the wind of the companion and dense 
mat erial following the co mpact object in its orbit 
(e.g. lHaberl fc Davlll992l) . The two bright >90 ks 
ObsIDs 6385, 6387 show an eclipse egress. Also 
ObsID 6377 and 6383 seem to indicate eclipses of 
PMH47. The frequency of the detection of eclipse 
in- and egresses in the 70 to 100 ks Chandra obser- 
vations points at an orbital period of the system 
that is not much longer than these observations. 

If we interpret the three transitions from low to 
high respective intermediate state in ObsIDs 6385, 
6387, and 6377 as eclipse egresses, the orbital pe- 
riod of the system can be determined. We use the 
fact that the times of egresses have to be sepa- 
rated by an integer number of orbital periods. We 
start from the shortest separation (ObsIDs 6385 
and 6377, separation r diff i = 6.931±0.026 d). Pos- 
sible orbital periods are P n = Tdiffi/n with n less 
than 7 (due to the shortest allowed orbital period 
of 1.04 days, see above). We now can test for which 
of these six candidate periods the time difference 
between the well determined eclipse egresses of 
ObsIDs 6387 and 6385 (T diff2 = 84.884 ± 0.018 
d) is consistent with an integer. Fortunately, we 
find a unique solution. For n = 4, Tdiff2 is - 
within the errors - 49 times P n and we get a 
best period from the ChASeM33 observations of 
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PchASoM33 = 1.73 2 33 ± 0.00037 d. 

This period can be further constrained by in- 
cluding XMM-Newton observations. Extrapolat- 
ing the eclipse egress from Chandra ObsIDs 6387 
with PchASoM33 to the time of XMM-Newton Ob- 
sIDs 0102641001 and 0102641101 (1096 periods 
earlier) shows that an eclipse egress should have 
happened between these observations (separated 
by 0.046 d) as the source is not detected until the 
end of ObsID 0102641001, but it is already bright 
at the beginning of ObsID 0102641101. Com- 
bining XMM-Newton and Chandra eclipse egress 
epochs we obtain an improved orbital period for 
PMH47 of P = 1.732479±0.000027d. The eclipse 
duration corresponds to 0.170 ± 0.007 in phase or 
an eclipse half angle of 30.6 ± 1.2 degree. 

In Tables Q] and [5] we give for the individual 
observations binary phase and cycle number at 
the beginning and end of the observations based 
on the best ephemeris determined above. It is 
clear that not all the variability can be explained 
by eclipses and additional variability has to be 
present. On the other hand some of the short 
XMM-Newton observations entirely fall in eclipse. 
Also the time without photons in Chandra ObsID 
6383 (see Sect. [4j coincides with an extrapolated 
eclipse. 

The ROSAT HRI observations are spread over 
12 days and cover several orbital periods of 
PMH47 (Fig. P. While during the first half of 
the observations PMH47 was detected most of the 
time, the rates are consistent with zero later on. 
A ROSAT HRI count rate of 1.0 x 10~ 3 ct s" 1 
corresponds to an un-absorbed luminosity in the 
0.2-4.5 keV band of 9xl0 36 erg s _1 assuming an 
absorbed power law spectrum as observed during 
Chandra ObsID 6387 (Ah = 6xl0 20 cm~ 2 , photon 
index T = 0.77). Therefore, during the ROSAT 
bright state, the PMH47 luminosity was ~ 5xl0 37 
erg s _1 , about a factor of two brighter than during 
the Chandra and XMM-Newton bright state de- 
tections. In Fig. [5]we mark the extrapolated times 
of eclipse using binary ephemeris derived from the 
Chandra and XMM-Newton data. The ROSAT 
count rates are plotted versus binary phase in 
Fig. [U The plot shows intensities compatible with 
zero at phases below 0.1 and above 0.85 indicat- 
ing that the extrapolation of the ephemeris by 
2560 periods still gives acceptable results. The 
extrapolation of the eclipse egress to the ROSAT 



observations leads to a phase uncertainty of 0.04. 
Unfortunately, the sampling of the phase space by 
the ROSAT observations during the bright time is 
not dense enough to determine an accurate eclipse 
egress which would allow us to improve on the bi- 
nary ephemeris. 

7. Optical counterpart 

The improved position of PMH47 (see Fig. 
II) suggests the star LGGS J013236. 92+303228.8 



dMassev et al.l 120061 ) to be the optical counter- 



part listed with 21.011 ± 0.011 mag in V with 
colors B - V of -0.122 ± 0.014 mag, U - B of 
-0.964±0.012, V-R of -0.062±0.018, and R-I 
of 0.033 ±0.014. 

This star was not identified as a variable in 



the M 33 CFHT variability survey (jHartman et al 
20061 ). To search for any small optical modulation, 



we derived the light curves of this star in the Sloan 
and i' bands from t he CFHT images usin 



.9 



(200< 



the procedure described in lHartman et al 
see Sect. [2]). Measured magnitudes are given in 
Table |2 

Folding with the X-ray period the 34, 33, and 
35 measurements in g', r', and i' bands, respec- 
tively, reveals a clear modulation. The folded light 
curves (see Fig. [7]) are of double sinusoidal shape, 
with one of the minima at the phase of mid X-ray 
eclipse, suggesting an ellipsoidal modulation of a 
high mass optical companion. The detection of an 
optical modulation at the X-ray period confirms 
the identification of the optical counterpart. 

We independently searched the optical data for 
periodicities. We constr ucted periodograms usin 
the method described at IShporer fc Mazehl (|200' 
Sect. 3.3) where for each frequency, each band 
was fitted with two harmonics and a zero point. 
The periodogram was taken to be the amplitudes 
sum, in quadrature, divided by the \ 2 . When us- 
ing all data points the 1.7325 day peak is the high- 
est, but there are many other similar strong peaks. 
When removing the data for the two outliers close 
to phase 0.35 and phase 0.53 (see Fig. [7]), there are 
only three significant peaks (0.0771 d , 0.5772 
d" 1 , 0.1517 d" 1 ), and the 1.7325 day peak is the 
second strongest (Fig. [5]). We checked if there is 
a justification for removing these points from the 
light curve based on sky conditions in the indi- 
vidual nights or image artifacts. While we can't 
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guarantee that there isn't some non-obvious sys- 
tematic variation, we are inclined to believe that 
the variations are real. 

8. Modeling the system 

In determining the parameters of the PMH47 
system, we first search for a stellar model (iden- 
tified by the effective temperature T e ff, the log- 
arithm of the surface gravity g, and the chemi- 
cal composition Fe/H in solar units) that fits the 
broad-band photometry of the system, assuming 
the optical brightness is coming from the opti- 
cal star alone. From the stellar radius, the ob- 
served X-ray eclipse width and the periodic ellip- 
tical modulation in the optical, we can derive the 
two masses as a function of the orbital inclina- 
tio n. This analys i s is si milar to the one presented 
bv lShporer etafl lj2007h for X-7. 



8.1. Optical companion 

To estimate the radius of the optical compan- 
ion we co mpared the UBVRI p hotometry for the 
star from iMassev et aD (j2006| ) to the [Fe/H] = 
table of model st ellar atmosphere br oad-band 
colors computed by iBessell et al.l (1998). Assum- 
ing a photometric uncertainty of 0.05 mag, color 
excess ratios of E(U - B)/E{B - V) = 0.71, 
E(V — R) I E(B - V) = 57 and E(R - I) /E(B - 
V) = 0.74 (|Bessell et al.lll998l) . and a minimum 
extinction of E(B — V) = 0.07 w hich is the fore- 
groun d reddening to M 33 from Ivan den Bergb 



(2000), we place a 95% lower limit on the effec 
tive temperature of 19000 K and a lower limit 
on log(<7) of 2.5. The formal best fit model has 
T cff = 33000 K, log(c/) = 4.5 and E{B-V) = 0.193 
with x 2 = 2.4. Note that the extinction is in good 
agreement with the range of values predicted from 
analyzing the X-ray spectrum (0.12 < E(B — V) < 
0.28). Assuming a distance of 795 kpc to M 33 
and the "standard value" for the ratio of total to 
selective extinction R v — Ay/E(B — V) = 3.1 
( Cardelli et al. 1989f ). we place a 95% upper limit 
of 10.8 Rq on the radius. For the best fit model 
the radius is 8.0 Rq. For the assumed distance 
and best-fit extinction the absolute magnitude of 
the star is My ~ —4.1 which, assuming the star 
is a main-sequence d warf, would correspo nd to 
a star of type B0V (|Schmidt-Kaleri Il982t ). We 
note that the effective temperature for such a star, 



T e f f ~ 30000 K, is close to the formal best-fit tem- 
perature from modelling the broad-band colors. 

8.2. Binary system parameters 

We used the " PHysi cs Of Eclipsing BinariEs 
(PHOEBE)" program of lPrsa fc Zwitteri 1)20051 ) . a 
front-end code for the Wi l son-Deyinney prog ram 
(jWilson fc Devinnevlll97ll : IWilsonlll979tll990h . to 



model the obtained periodic light curves and to 
derive an estimate of the masses of both compo- 
nents. The main effect is due to the tidally in- 
duced ellipsoidal shape of the optical component. 
We ran PHOEBE on the CFHT g',r', and i' light 
curves, assuming the X-ray determined eclipse half 
angle value of 30.6° ± 1.2 ° and using l inear limb 
darkening coefficients from lClaret ( 2004 ). stepping 
through inclination and fitting for the mass ratio. 
Figure [7] shows a representative PHOEBE model 
fit to the light curves for a binary system with 
an inclination of i — 80° . We note that due to the 
degeneracy between q, i and the Roche-lobe-filling 
factor, the model fit is similar for other allowed in- 
clinations. Figure [9] shows a plot of the optical star 
and compact-object masses as a function of i for 
each of these different T e g values. Below i = 70° 
the optical star fills its Roche Lobe. From the plot 
for the compact object mass, a neutron star with 
canonical mass of 1.4 M© seems to be allowed for 
inclinations below 85° for the 50000 K model, but 
only for below 75° for the 33000 K model. For 
the 19000 K model, a canonical mass neutron star 
seems to be excluded. A neutron star with canon- 
ical mass at the best fit companion temperature 
T e ff of 33000 K would imply a system inclination 
i = 72° and a companion mass of 10.9 M@. 

9. Discussion and summary 

The time variability of PMH47 in the ChASeM33 
observations can be explained with a 1.7 d eclips- 
ing X-ray binary with an eclipse half angle of 
30.6°. The identification with an OB star in M 33 
that shows an ellipsoidal light curve, clearly estab- 
lishes the system as the second HMXB in M 33. 
We do not detect pulsations in the X-ray flux and 
therefore the compact object in the system could 
be a neutron star or a black hole. The hard X-ray 
spectrum, howe ver, may indicate a neutron star 



system (see e.g. lWhite et al.lll983l ). 

We can not exclude that the long term variabil- 



es 



ity of the system in the Chandra, XMM-Newton 
and ROSAT observations can be of irregular na- 
ture. However, it also could reflect a long term pe- 
riodicity similar to those known in Her X— 1, LMC 
X-4 or SMC X-l. The Chandra bright state ob- 
servations (ObsID 6387 and 6385) are separated 
by 84 days. With our sparse unsystematic moni- 
toring, we can not decide if the source had another 
on-state in between as would be expected by ex- 
trapolating from the 35 d period of Her X-l, a 
system with a similar orbital period. 

1 2006h 



Inspecting the available iMassev et al 



colour and emission-line images we find that the 
object is located in a region with few (bright) 
stars. Also, there are no close-by emission nebulae, 
indicating that the object is older than M 33 X— 7 
which is still located in a H n region. 

PMH47 was serendipitously detected far off axis 
in the ChASeM33 observations. Only with the 
help of dedicated deep observations with good 
time resolution it will be possible to characterize 
long term variability of PMH47 and to search for 
the expected pulsation period below 50 seconds. 
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Table 1: Chandra observations of the ChASeM33 program covering [PMH2004] 47. 



ChAScM33 Obs.id. Obs.dates a On-time ACIS Offax [PMH2004] 47 binary 



Field (ks) CCD-ID (') phase b cycle b comment 





1730 


2000-07-12 


40 c 


7 


17.6 


0.795-1.062 


-1305,-1304 


ingress, off 


4 


6382 


2005-11-23 


73 





10.2 


0.168-0.589 


-173 


faint 


4 


7226 


2005-11-26 


25 





10.2 


0.703-0.856 


-172 


faint 


4 


6383 


2006-06-15 


99 


3 


10.8 


0.754-1.424 


-56,-55 


eclipse, bright 


6 


6387 


2006-06-26 


78 


6 


19.7 


0.993-1.507 


-50,-49 


egress, bright 


6 


7344 


2006-07-01 


22 


6 


19.7 


0.882-0.998 


-47 


ingress, off 


5 


6385 


2006-09-18 


91 


2 


10.9 


0.534-1.139 


-1,0 


eclipse, bright 


1 


6377 


2006-09-25 


94 


6 


17.4 


0.906-1.507 


3,4 


egress, faint 



"date of start of observation 

6 with respect to eclipse center HJD 245 3997.476 and orbital period 1.732479 d (see text) 
c rcsidual exposure after screening first 10 ks due to background flaring 



Table 2: XMM-Newton observations covering [PMH2004] 47. 



Obs.id. 


Obs.datcs a 


On-time 


Off ax 




[PMH2004] 47 binary 








(ks) 


(') 




phase c 


cycle c 


comment 


0102640401 


2000-08-02 


16 


5.6 


< 33 


0.952-1.039 


-1293,-1292 


eclipse 


0102640501 


2001-07-05 


12 


10.7 


836 ± 65 


0.729-0.809 


-1098 


edge of pn, bright 


0102641001 


2001-07-08 


12 


14.7 


< 206 


0.992-1.072 


-1097,-1096 


eclipse 


0102641101 


2001-07-08 


12 


3.2 


1732 ± 68 


0.099-0.179 


-1096 


bright 


0102642301 


2002-01-27 


13 


13.7 


< 67 


0.305-0.385 


-979 


off 


0141980601 


2003-01-23 


11 


5.8 


< 13 


0.919-1.011 


-771,-770 


eclipse 


0141980101 


2003-07-11 


16 


14.7 


< 205 


0.377-0.484 


-673 


off 



"date of start of observation 

'0.2-4.5 keV absorption corrected luminosity or 3<r upper limits in units of 10 34 erg s —1 assuming an absorb ed power law 
spectrum (Afjj = 6xl0 20 cm -2 , photon index T = 1.7) and a distance to M 33 of 795 kpc llvan den Bergblll99iri which we use 
throughout the paper 

c with respect to eclipse center HJD 245 3997.476 and orbital period 1.732479 d (see text) 

Table 3: Power law spectral modeling results for [PMH2004] 47 for the bright time in the Chandra ACIS- 
I ObsIDs 6385 and 6387 and for the XMM-Newton EPIC observation 0102641101 (all EPIC instruments 
fitted together). For each instrument, we give the effective integration time ti nt and the raw count rate. The 
number of energy bins reduced by the number of free parameters defines the degrees of freedom v. 90% 
errors are given. 



Observation 


Inst. 


tint 

(ks) 


Rate a 


-^HM33 b 


r 


L x c 


// 




Chandra 6385 


ACIS I 


63.1 


2.27 


iD -12 


0.94 ±0.13 


1.9 


66 


1.38 


Chandra 6387 


ACIS I 


62.5 


2.20 


< 8 


77+° 12 
u -' '-0.10 


2.0 


77 


1.19 


XMM-Newton 
EPIC 

0102542301 


PN 

MOSl 

MOS2 


8.8 
9.8 
9.8 


6.40 
1.73 
1.77 


-i 9+8.5 


0.88 ±0.13 


1.8 


28 


1.38 



"Raw count rate in units of 10 -2 ct s _1 as given in XSPEC 

6 Absorption in units of 10 20 cm -2 exceeding the fixed Galactic foreground of 6.0xl0 20 cm -2 

c X-ray luminosity in the 0.2-4.5 keV band in units of 10 37 erg s — 1 , corrected for absorption, for extraction radii, and for 
vignetting 
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Table 4: Optical light curve of the star LGGS J013236. 92+303228.8 identified with [PMH2004] 47. Magnitude 
with error in the g' , r' , and i' filter and heliocentric Julian date (HJD) of the observation are given. 



HJD 






g' 




HJD 






r' 




HJD 






i' 


245 0000.0 


(mag) 


245 0000.0 


(maj 





245 0000.0 


(mag) 


2873 


984375 


21 


077 


'13) 


2873 


992665 


21 


424 


;i7) 


2874 


001810 


21 


731(39) 


2876 


041448 


21 


006 


'12) 


2876 


049696 


21 


324 


16) 


2876 


058663 


21 


606(35) 


2876 


067077 


20 


982 


'12) 


2876 


075340 


21 


308 


17) 


2876 


086306 


21 


630(36) 


2881 


949358 


21 


039 


'15) 


2881 


957628 


21 


346 


17) 


2881 


966625 


21 


668(36) 


2882 


917645 


20 


999 




2882 


925921 


21 


321 


19) 


2882 


934887 


21 


673(36) 


2882 


943189 


20 


987 


'12) 


2882 


951458 


21 


317 


16) 


2882 


960463 


21 


645(34) 


2886 


083231 


21 


085 


'14) 


2886 


091515 


21 


430 


18) 


2886 


100505 


21 


720(40) 


2903 


847418 


20 


973 


'14) 


2903 


855644 


21 


328 


18) 


2903 


864595 


21 


604(36) 


2905 


975845 


21 


050 


'14) 


2905 


984073 


21 


374 


19) 


2905 


993022 


21 


707(43) 


2908 


917034 


20 


982 


'13) 


2908 


925269 


21 


334 


16) 


2908 


934226 


21 


616(32) 


2910 


026011 


20 


988 


'15) 


2910 


034266 


21 


330 


19) 


2910 


043261 


21 


622(51) 


2931 


777952 


21 


031 


'14) 


2931 


786238 


21 


376 


'20) 


2931 


795225 


21 


560(42) 


2931 


804198 


21 


035 


'14) 


2931 


812455 


21 


369 


'20) 


2931 


821460 


21 


703(43) 


2931 


829883 


21 


055 


'12) 


2931 


838149 


21 


407 


18) 


2931 


912971 


21 


720(44) 


2931 


921486 


21 


046 


:w) 


2931 


929751 


21 


389 


19) 


2931 


938758 


21 


708(38) 


2940 


741715 


21 


026 


;i4) 


2940 


749948 


21 


376 


'22) 


2940 


758913 


21 


713(51) 


2963 


800785 


20 


940 


;i4) 


2963 


809086 


21 


267 


'21) 


2963 


818136 


21 


477(41) 


2991 


813734 


20 


981 


;i4) 


2991 


821989 


21 


302 


'20) 


2991 


830949 


21 


530(40) 


3240 


978670 


20 


988 


;is) 


3240 


986361 


21 


335 


17) 


3240 


994707 


21 


616(37) 


3256 


073787 


21 


042 


;i6) 


3256 


081416 
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Fig. 1.— left: Chandra ACIS-I light curves (0.5-5 keV) of ObsID 1730, 6382, 7226, 6383, 6387, 7344, 6385 
and 6377 integrated over 5000 s, 5000 s, 5000 s, 5000 s, 1000 s, 5000 s, 1000 s and 3000 s, respectively. Time 
zero corresponds to HJD 245 0000.0 + (1737.8527, 3698.0480, 3700.7070, 3901.7629, 3912.5724, 3917.5767, 
3996.6692q and 4004.2427), respectively (solar system barycentre corrected), right: 6' x 6' ACIS images 
(0.5-5 keV) of the [PMH2004] 47 area for the individual ObsIDs. Source and background extraction regions 
are indicated by rectangles and ellipses. The dashed ellipses show the HRMA PSF 90% enclosed count 
fraction regions. 
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Fig. 1. — Continued. 
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Fig. 3. — Chandra ACIS-I spectrum of 
[PMH2004] 47 for the bright state of ObsID 6387. 
The histogram shows the best-fit absorbed power 
law model. 




o 30 - 



70 75 80 85 90 
Inclination [Degrees] 




S 



70 75 80 85 90 
Inclination [Degrees] 




Fig. 9. — Plots showing the mass vs. inclination 
for three values of the stellar temperature. The 
left plot is for the optical component, the right 
plot is for the compact object. The solid, the 
dotted and the dashed lines show the results for 
T cff =33000, 19000, and 50000 K, respectively, us- 
ing the PHOEBE code. 



Fig. 8. — Multi-band double-harmonic peri- 
odogram after removing outliers around phase 
0.35 and 0.55. Peak 2 corresponds to the X-ray 
period (see text). 
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Fig. 2. — 6' x 6' Chandra ACIS-I 0.5-5 keV images of observation 6385 for times before the ingress of 
[PMH2004] 47 into low state (left), during low state (middle), and after low state (right). The position of 
[PMH2004] 47 is indicated by the ellipse. 
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Fig. 4. — Local Group survey B image (|Massev et al.ll2006f ) of the [PMH2004] 47 field. X -ray positions are 
i ndicated by circles wit h radii of the total 3er error. The circle to the SW gives the XMM- Newton position 
I Misa novic et al. I l2006h . The dashed circle to the NE indicates the ChASeM33 position (jPlucinskv et al 



2008). The small circle indicates the improved C hASeM33 position ( see text). The bright object included 



in all error circles is LGGS J013236.92+303228.8 (|Massev et alj |2006). 
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Fig. 5.— ROSAT HRI (0.1-2.4 keV) light curve of [PMH2004] 47 from July 27 to August 8, 1994. Times of 
eclipses are indicated using eclipse center HJD 245 3997.476 and orbital period 1.732479 d. 
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Fig. 6.— ROSAT HRI (0.1-2.4 keV) count rate of [PMH2004] 47 versus binary phase (eclipse center 
HJD 245 3997.476 and orbital period 1.732479 d, see text). Phase of eclipse is indicated. 
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Fig. 7. — Folded light curves, in magnitude, in (top to bottom) g' , r' and i' . All light curves are folded using 
the X-ray ephemeris. The dotted line shows a PHOEBE fit for inclination 80° and effective temperature of 
33000 K. 
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